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^ I Abstract. We shall show that the density of defects produced at a second- 

' order phase transition is determined by the correlation length of the fields. 

I ' This is true both for defects appearing in the Ginzburg regime and for 

2 ■ defects produced at a quench, when the Ginzburg regime is irrelevant. 

a 



^ . 1. INTRODUCTION 

These notes are based on lectures given by one of us (R.J.R) at the Nato Advanced 
Study Institute and Euroconference on Formation and Interactions of Topological De- 
fects held at the Isaac Newton Institute, Cambridge (UK) in September, 1994. They 
aim to provide a preliminary discussion of how strings (vortices) and monopoles can be 
produced at the phase transitions of relativistic quantum fields. 

The applications that we have in mind are to the early universe, where it has been 
argued (see Kibble, these proceedings and elsewhere [|l|) that cosmic strings produced 
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at the era of Grand Unification (i.e. at energy scales 10^^ - 10^^ GeV) can provide the 
seeds for the large-scale structure formation in the universe that we see today. There 
are many attractive features to this idea, and we shall not recount them here. However, 
it should be remembered that the isotropy of the universe suggests that it has passed 
through a period of rapid inflation. Inflation in itself generates large-scale structure 
but, of greater importance in this context, it was originally introduced to dilute the 
undesirably high monopole density that is almost inevitable in unified theories. Some 
ingenuity is required for strings not to suffer a similar fate. Should cosmic strings turn 
out to be nothing more than an elaborate fancy we are consoled by the observation that 
they are, in many ways, the relativistic counterparts to the vortices in superfiuids and 
superconductors. So much so, that vortex production |^ in superfiuid "^He has been 
invoked as simulating cosmology in the laboratory (see Zurek, these proceedings, and 
elsewhere ^ ). But for their relativistic nature, our methods are equally applicable to 
these more homely materials. [By homely, we mean terrestrial rather than non-exotic. 
See Salomaa, these proceedings, for a description of some of the very exotic defects of 
superfiuid ^He, for example]. However, we shall not change to a non-relativistic gear 
here. 

In practice, our calculations yet, too primitive to be able to address the 

details of early universe cosmology directly, even if we had wished. However, our main 
conclusion is general, and does not require a cosmological backcloth. It is that defects 
play an important role in any second-order transition at which they can be produced, 
appearing initially in essentially the maximum numbers compatible with retaining their 
individuality as diffuse entities. Specifically, the density of defects, after their produc- 
tion at the phase transition, is determined largely by the distance ^{t) over which the 
fields are correlated, and which also characterises defect size. Very crudely, we predict 
one monopole per correlation volume v = 0{C,^{t)) and one string passing through each 
correlation area a = 0{^'^{t)). This applies equally to 

i. defects produced from large fluctuations in the Ginzburg regime near a second- 
order transition and to 

ii. defects produced from long wavelength spinodal decomposition after a quench. 

In the first case this result was anticipated qualitatively by Kibble, on making rea- 
sonable assumptions about domain formation. It has now become sufficiently part of 
the folklore that the results in the second case will be no surprise, even though circum- 
stances are different there. However, our conclusions are reached as a result of quantum 
calculations rather than semiclassical arguments which, for gauge theories, have led to 
some confusion. Of course, there are several caveats (e.g. weak couplings, short to 
intermediate times) as will be seen. Whether these would be satisfied in realistic many- 
body systems is dubious, but our calculations can be improved upon in principle. In 
practice we have yet to do so, and the results presented here are chosen, in part, for 
their simple analytic nature. 

These lecture notes essentially fall into two parts, corresponding to cases i) and 
ii) above. The first is concerned with phase transitions as viewed from the platform 
of equilibrium thermal field theory. To be concrete and simple we restrict ourselves 
to the transitions of global and local 0{N) theories. In D = 3 spatial dimensions an 
0{N) theory permits strings when N = 2 and monopoles when = 3. We examine 
the presence of both global and local defects in the Ginzburg regime close to a second- 
order transition, characterised by large fluctuations. However, while confirming our 



basic ideas about fluctuations we are unable to provide a satisfactory mechanism as 
to how defects appearing in this way can persist as the system cools. The second 
part of the lectures attempts to circumvent this problem by avoiding any discussion 
of equilibrium theory, adopting a non-equilibrium approach from the start for global 
0{N) theories. Specifically, we consider the production of defects as a consequence 
of spinodal decomposition at the onset of a transition, whose freezing in is much less 
problematical. Even then, a quantitative description has yet to be given. We conclude 
with some tentative steps in this direction. 

Our introduction to equilibrium and non-equilibrium quantum field theory is, or 
should be , well-known. Some of the further material is taken from unpublished lecture 
notes (T.S.E) and from published work (R.J.R) with Mark Hindmarsh |Q. The work 
on defect production is more recent, as yet unpublished (although in a preliminary 
form some of the conclusions were presented P] (R.J.R) at the recent Nato Advanced 
Workshop on Electroweak Theory and the Early Universe, held in Sintra, Portugal, 
March, 1994). However, the main results of the latter part have already been submitted 
(R.J.R. and A. Gill) §. 

In addition we have drawn heavily on the recent series of papers by Dan Boy- 
anovsky. Hector de Vega and co-authors [0, |^, which provide an excellent introduction 
to behaviour out of equilibrium, and to the much earlier (equilibrium) work of Halperin 
concerning global defects. We are indebted to Mark Hindmarsh for helpful comments 
on the latter. The reader looking for background material to amplify some of our more 
cursory comments should find the Proceedings of the Third Thermal Fields Workshop, 
Banff (Canada), August 1993, most helpful [TO. 



2. THE PHASE TRANSITIONS OF RELATIVISTIC QUANTUM 
FIELDS 

The defects in which we are interested are not fundamental entities like superstrings. 
They are diffuse field configurations formed at phase transitions in the early universe, 
which survive because of their topological stability. To understand how phase transi- 
tions occur we need to recapulate the rudiments of thermal field theory |jlT|, |12|. 

2.1. What is Thermal Field Theory? 

Thermal field theory is a combination of two theories. The first is relativistic quantum 
field theory, used to describe the behaviour of elementary particles, when only a few 
such particles are involved e.g. in — e~ collisions at LEP. The specific attribute that 
special relativity brings is the annihilation and creation of particles, the conversion of 
rest-mass to energy. This is combined with the fundamental ingredient of quantum 
physics, quantum fluctuations as encoded by Heisenberg's uncertainty relation 

AEAt > h. (1) 

One of the great steps forward in the development of quantum field theory was 
the use of Feynman diagrams to represent the effects of quantum fluctuations on phys- 
ical processes. For instance, for the case of an electron travelling in vacuo, in quantum 
electrodynamics we must include corrections such as the first diagram in Fig. |l] The in- 
ternal lines represent interactions between the electron (the external legs) and possible 




Figure 1. Different types of process which occur when in a heat bath. 



vacuum ffuctuations involving the emission and absorption of photons (the wavy hne). 
The latter are virtual, as are the intermediate electrons, existing only for the short 
time allowed by Heisenberg's uncertainty principle. The picture we have is of an elec- 
tron travelling through a vacuum, empty of real particles but full of virtual shortlived 
particles, which interact with it. 

The second ingredient in thermal field theory is thermodynamics/ statistical me- 
chanics, used to study many-body problems. The key idea is that a few thermodynamic 
or bulk properties are sufficient to characterise the essential physics. In doing this, the 
precise initial state of the system is assumed unknown, the choice of any particular 
state only specified by a probability. This statistical uncertainty is quite different from 
the quantum uncertainty discussed above. In particular, we shall typically assume a 
Boltzmann distribution, e~^^, to describe the probability of being in a state of energy 
E. In units in which Boltzmann's constant ks = 1, (3 is the inverse temperature 
of the system. When combined with quantum physics this leads to fluctuations in the 
number of real particles, as described by the usual Bose- Einstein or Fermi-Dirac number 
distributions ^ 

= (2) 

(with =p taken according as the particles are bosons/fermions). A useful picture is to 
think of doing experiments in a background, the heat bath or reservoir, which is full of 
real particles, whose precise number is unknown. 

This combination of relativistic quantum field theory and statistical mechanics, 
which constitutes thermal field theory, describes quite a different physical situation from 
that of normal quantum field theory. Both quantum and statistical fluctuations have 
to be accounted for simultaneously and, as outlined above, they are dissimilar in their 
nature. Although it is not surprising that thermal field theory has some quite different 
properties from standard quantum field theory, what is amazing is that the two can be 
described in a very similar way. 

As an example of thermal physics, return to the electron, now taken to be prop- 
agating through a QED plasma. The quantum fluctuations mentioned above (the first 
diagram in Fig. 0) are always present, whereas the remaining diagrams in the Figure 
do not appear in ordinary quantum field theory, involving interactions with the real 
particles present in the plasma or heat bath. The last diagram in Fig.Q represents such 
an interaction with two real particles, while the middle two diagrams represent the 
effects on the electron's propagation of a mixture of quantum (virtual) and statistical 
fluctuations. All of these will change the electron's inertial mass, and the latter lead to 
the dissipation of its energy. 



2.2. Why Are There Phase Transitions? 



The diagrams of Fig. |T] suggest that a particle in a heat-bath behaves as if its 



momentum-space propagator were [ITS 



Gik) = + 2n6{e - w?)n,{u). (3) 

The first term is the usual virtual particle exchange, present for all four-momenta fc, 
whereas the second term describes real particles in the heatbath, present only when 
k"^ = m? (in units in which c = 1). We have taken spinless bosons as an example. It 
will be shown later that the situation is rather more complicated, but G{k) of (2.3) is 
all that is necessary for one-loop diagrams, which are sufficient to show the presence of 
phase transitions. 

Consider the theory of a real scalar field 0, given by the action 

S[<P\ = J d'x[^-{d,<j>){dy) - \ml<p' - \gl<p'] (4) 

in thermal equilibrium at temperature T. The presence of the heatbath affects the 
inertial properties of particles fired into it. The effective mass m^(T) of (j) -field quanta 
is represented diagramatically (to one loop) by 





+ (5) 



where the solid line now denotes the full propagator (3), and the first term is mg, read 
off from (4). [The fact that m^(T) is independent of momentum is a refiection of the 
relative simplicity of the one-loop diagrams in g'^ip'^ theory, in contrast to those of QED.] 
The second term of (5) is ultraviolet divergent, but this causes no problem, since the 
bare mass mo has no physical meaning. To O(^), the physical (mass)^, denoted by m?, 
is defined as the sum of the first two terms. The end result is that the heat-bath induces 
a temperature-dependent term to the effective mass of the form 

3 
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where cf^/c = {2Tc)^^d^k, and Lu(k) = yfc^ H-m^. Since the Boltzmann factor cuts off 
the momentum integration at \k\ = 0{T) the second term in (6) is ultraviolet finite. Its 
temperature-dependence (for T ^ m) can be read off from its ultraviolet divergence, 
had nfe been absent, as 0{g'^T'^). Specifically, up to logarithms in the final term 

m\T) =m^ + -hg\T^ --Tm + 0(m^)). (7) 
8 vr 

[At one loop there is no coupling constant renormalisation and we have renamed go as 
g]. That is, the effect of the heatbath is to increase the (mass)^ value of the quanta. 
In particular, massless bosons put in a heatbath acquire mass. With qualifications this 
is true for photons, as will be seen later. 



Suppose now that the cold (mass)^ parameter is negative, rn? = —\g'^rf , say, 
some T] . That is, on redefining the parameters to accomodate the one-loop quantum 
fiuctuations, the action S'[0] for is 

S[<P\ = J d'x[l{d,<j>W<j>) - - (8) 

the familiar double-well potential, with vacuum values (p = ■ The choice of one 
of these minima as the ground state breaks the refiection invariance (p ^ ~4' of 
action (or Hamiltonian) and this Z2 symmetry is said to be spontaneously broken. 
The upturned quadratic potential —\g^Ti'^(p'^ near the origin characterises a region of 
instability for field fiuctuations of small amplitude. For high T, when only the first 
term in the brackets of (1.8) need be taken, m?{T) becomes positive, however negative 
may be. That is, the Z2 symmetry is restored |T^, |1^]. The effective (mass)^ 



for small field amplitudes is then most transparently written as 

m\T)=m'[l-^)j, (9) 

where T'^ = At]^ and we are working in units in which h = 1. In this approximation Tc 
defines the critical temperature of a second-order transition, above which the symmetry 
is restored, and below which it is broken. At the transition the correlation length 
^(T) = \m(T)\'^ diverges as 

e(r)(x|r-r,|-^7 = i (lo) 

This mean-field result can be improved by a renormalisation-group analysis but, at the 
qualitative level at which we are working, it is good enough. 

2.3. Phase Transitions of 0(N) Theories. 

The theories in which we shall be interested are particular examples of the 0{N) ex- 
tension of the theory (4), in which N scalar fields (pa (« = 1,2, ...A^) transform as the 
fundamental vector representation. If we are only concerned with invariance under 
global 0{N) transformations, the action (4) is elevated to 

S[<p] = I d'x[lid,<PaW<Pa) - Ig'i^Pl - v'n (11) 

where D = 3 is the number of spatial dimensions and summation over the 0{N) label a 
is assumed. For rj^ > the 0{N) symmetry is broken to 0(A^ — 1), the vacuum manifold 
A4 = 0{N)/0{N — 1) being equal to S*^^^, the (A^ — l)-sphere. On expanding about 
any point on the sphere we find a single Higgs boson with mass ttlh = gf] and (A^ — 1) 
massless Goldstone bosons. If we heat the system to a plasma at temperature T we 
would expect that, as for the single field, at some temperature = 0{7]) the symmetry 
would be restored. This is indeed the case. Some care is needed to guarantee the 



masslessness of the Goldstone modes at all T < Tc [|12| but basically we proceed as 
in (6) in determining the effective mass m{T) of the small-field potential. The only 
difference is that the loop diagram of (5) is extended to all N scalar fields. The end 
result is that, ignoring terms relatively 0{r]/Tc), the one-loop calculation gives |ll3| 



= —gW + (N + 2)c/^ (12) 



with = 12rf/[N + 2). For T < Tc the symmetry is broken, with global minima at 
101 = viT), where r]^{T) — r]'^{l — T'^/T^). The effective Higgs mass mn, that measures 
the curvature in the radial field at these minima, is mniT) = gT]{T) or,equivalently, is 
given by rn\j{T) = — 2m^(T). As before, there is a second-order transition (with the 
same index 7 = ^ at one loop). 

Global invariance sits uneasily in contemporary particle physics, for which there is 
no evidence for Goldstone bosons. To extend the theory to be invariant under local 0{N) 
transformations, it is necessary to introduce \N{N — 1) gauge fields transforming 
as the adjoint representation of 0{N). The Lagrangian density can be written as 

C = \D^^D^ct> - -TrF.^F^'' - 1/(0^ - r,^. (14) 
Here the covariant derivative of is 

D^,(t> = + eA^(t>, (15) 

and 

F^. = d^A, - d,A^ + e[A^, A,]. (16) 

The {N — 1) Goldstone modes of the global theory now transmute into longitudinal 
modes of the vector fields, enabling [N — 1) gauge fields to acquire mass = er], while 
the Higgs mass is mn — gv^ before. The remaining gauge fields stay massless. In 
calculating ■m^{T) we now have to include gauge-field (wavy line) one-loop diagrams 





The precise diagrams involved depend on the gauge (these are appropriate for the 
Landau limit of the covariant gauges C = 0) some of these diagrams do not give 
O(T^) contributions. Further, we have not included tadpole corrections so that we are 
assuming that we have shifted the vacimm expectation value 1] vi'^)- At leading 
order inT/rj , the {N — 1) massive gauge modes cause the effective w?{T) of the global 
scalar theory to be changed from (12) to 



[6(iV - lje' + (iV +2j^''^ 
lowering the temperature at which m?{T) vanishes to 



m\T) = -\gW + m - l)e' + (A^ + 2)/]^, (17) 



rp2 ^ ^'^aW ..ox 

" &{N -l)e^ + {N + 2)g^' ^ ' 

However, a second-order transition at temperature is no longer guaranteed. The 
approximation of retaining only the O(T^) term in (17) may not be valid if the ratio 
e/ g is sufficiently large. What we have calculated here is the \im?'{T)(j)l^ contribution to 
the scalar sector effective potential. Had we retained terms of relative order rriv/T, the 
one-loop contribution of the gauge field is not so much {N — l)e^T^/A, as seems from 
(17). More accurately, it is 

m\T) = 6{N - l)e' - ^Tmy(0)] , (19) 



where mv{(f)) = e|0| becomes mv when |0| = 77. This term, hnear in |0|, induces a 
cubic term O(e^T|0|^) in the effective potential |TB . 

Such terms, if strong enough (i.e. e/g large enough) can turn the second-order 
transition into a first-order transition at a temperature that is still essentially Tc of (18). 
Some caution is necessary. Eq. (19) is too simplistic. As will be seen later, the effect of 
the plasma is to induce electric and (for N > 1) magnetic screening masses to the gauge 
fields, and my (no longer a single mass) is not simply proportional to |0| |T6|. Further, 
magnetic screening is intrinsically non-perturbative in e. Nonetheless, for e/ g ^ 1, we 
do expect a first-order transition. Defect formation at a strong first-order transition 
is different from that at a second-order transition, proceeding largely by bubble nucle- 
ation. We shall not consider such a possibility and only assume second-order transitions 
(inevitable for the global theory, as we have seen). In fact, this may not be a great 
restriction. The satisfactory evolution of the early universe string network requires a 
mechanism for strings to chop one another up and for the fragment loops to decay. 
A strong first-order theory would have strings with very different intercommutativity 
properties from second-order strings (small e/g), and it could well be that they would 
have undesirable consequences. 



3. FIELD FLUCTUATION PROBABILITIES IN THERMAL EQUILIB- 
RIUM 

Topological defects like strings and monopoles can be characterised by non-local field 
configurations (e.g. magnetic flux through a surface) or by local configurations (e.g. 
field zeros). In either case, the likelihood of their appearance at phase transitions can 
be determined if we know the probabilities for arbitrary field configurations. 

We begin by calculating the configuration probabilities of 0{N) scalar fields in 
thermal equilibrium at temperature T. To do so it is more convenient to adopt the 
imaginary-time approach to thermal field theory than the real-time approach (with its 
heatbath populated by real particles whose propagation we follow) adopted earlier. 

3.1. The Imaginary-Time Formalism. 



The idea is simple [|TT|, |T^. First consider the theory of the single scalar field (j) of (4). 
Let H[7T,(j)] be the Hamiltonian derived from S[(j)] (with vr = <p). Then, at temperature 
T = j3~^, the partition function Z 

Z = trp = tre-^^ (20) 

can be written as 

Z = J2<^n,to\e-^^\^n,to> (21) 

n 

where, in evaluating trp we have chosen a basis of eigenstates of (f) at time to- i-e. 

(f){to,x)\<t>n,to >= <^n{x)\<^n,to > . (22) 

For simplicity, the $„ have been taken to be denumerable (e.g. by the imposition of 
periodic boundaries). 



Already, the interplay between thermal and quantum fluctuations is apparent. 
The diagonal matrix element < to|e~''^^|$Tn > permits two interpretations:- 

<^n,to\e-^"\^n,to>^Pto[^n], (23) 

the (relative) probability that the field takes the value $n(^) at time to, or 

< to\e~^"\^n: to > = < to - iPM^n: to >, (24) 

the probability amplitude that the field has value at time to — ifi if it had value 
at time to- We continue to work in units in which h — 1. 

Since the system is in thermal equilibrium both interpretations are independent 
of to- The first (statistical mechanical) is just a reiteration of the probabilities of the 
Boltzmann distribution and gives us the quantity p[<l>] we wish to calculate. However, 
our ability to interpret this probability as a probability amplitude permits the quantum 
mechanical path integral realisation 



J4>{to,x, 



V(f) e 



(25) 



where the sum is restricted to fields periodic in imaginary time, (p{to, x) = (f){to — iP, x), 
period /3. 

The periodicity of the field in Euclidean (imaginary) time is most simply imple- 
mented by integrating in a straight line from to to to — i/3 as in Fig. 2. 



h 



Figure 2. ITF curve. 

of (25) then becomes 

ptom 



V(f) e 



where 



SM - [dr I d'xil^' + 1(V0)^ + ^m^0^ + 



(26) 
(27) 



for fields periodic in r with period P (position labels have been omitted, and the dot 
denotes differentiation with respect to r). This periodicity is enforced by Fourier ex- 
panding 0(r, x) as 

0(r,f)=5:0„(f)e2™-/^ (28) 

n 

whereby the quadratic part of Se[(P] becomes 

Se[<I>],u = /5E / ^'^[KWn)' + >M (29) 

in which ^ 

9 9 / 27m \ , , 

< = ml+[-^j ■ (30) 

At high T the ratio mn/frio is large. The n = mode is the light mode of the field, the 
n 7^ modes the heavy modes. One very successful approach to high -T equilibrium 
theory is to integrate out the heavy modes, to obtain an effective t/iree-dimensional 0o 
- theory. That is, we rely on the controlled decoupling of massive modes to perform a 
dimensional reduction from four to three dimensions, valid when 

(31) 

where L is the system size, ^ the field correlation length, and (3 the thermal wavelength. 
In estimating probabilities, it is good enough to observe that, from (25), 

PtoM = /n^</'n^(E</'n(0)-$)e-^-[^l (32) 

n n 

~ J V<f)o 5(0o(O) - $) /n ^</>ne-^-['^l, (33) 



since heavy mode fluctuations are small. The 6 denotes a 5-functional in which 
0o(O,x) = $(af) at each point in space. On defining the three-dimensional effective 
action 6*3 [0o] by 

N^-PSsiM = fY[v<f)n e-^^[<^] (34) 

where iV is a normalisation constant, it then follows from (33) that 

PtJ*] ~ A^e-^^3l*°l (35) 

[A better calculation |^ would show that 

PiJ<l>] = Are-^^[*°l, (36) 

where 

/?2 ff5S-i'' 



^[*] = ^3[$] - ^ y (^^j +0{P') (37) 

but, at the level to which we are working, (35) is sufficient.] For S'£;[0] of (29), the 
one-loop high-T behaviour of 6*3 [$] (in which only terms 0{g'^T'^) are displayed) is 

= J d'x[l{V<^y + im2(T)$2 + (38) 

where m?{T) is the effective mass (9). It is now understood as the light-mode mass 
whose vanishing signals the occurrence of a phase transition. [By not containing 0o 
fiuctuations,m^(T) of (38) differs from m?{T) of (9) by non-leading terms, which we 
ignore]. A more accurate calculation of would show it to be non-local. The 

expression (38) is sensible only for distances larger than the thermal wave length /3, at 
which there is an effective cutoff. 



3.2. Gaussian Fluctuations. 

Remaining witli tlie single field, the probability Pfo[*^'], that (f)(tQ,x) = is just a 

starting point. From it we wish to calculate the more general probabilities p{Q[4>] = q) 
that some functional of 0(to; Qi^], equals q. The quantity q could be a number, or a 
function q{x) , depending on Q. For simplicity, take Q as a pure functional, q a number. 
Generahsation is straightforward. If < ... > denotes averaging with respect to Pto[^], 
then 

p{Q['i>] = q) = <S{Q[cf>]-q)> (39) 
= jda< e^^^^t-^l-^) > . (40) 

The cumulant expansion of (40) is, up to quadratic terms 

p(Q[0] =q)= Jda e--ye-^"'[<'2'2>-<^>'l+-. (41) 

In the Gaussian approximation, in which only these terms are retained, the a integration 
can be performed, to give 

p(g[0] =q)= Are-5(9-<Q»'/<W>c, (42) 

where < QQ >c=< QQ > — < Q >^ is the connected two-Q correlation function. 

Our immediate interest is to determine in what circumstances field fiuctuations 
are large. The simplest measure of fluctuations is the coarse-grained fleld Ql'Cp] = (pv '■- 

0^ = -/ d^x(P{to,x) (43) 
V Jxev 

for an arbitrary volume v. [We have used the symbol v to denote both the position of 
the volume and its magnitude]. Introducing the window function (indicator function) 
I{x): I{x) = l,x G v; I{x) = 0,x ^ v, it follows that 

< (f)y(f)y >= \ [ dxdx'I(x)I(x')G(x - x'), (44) 

where G{x — x') is the equal-time correlation function, read off from (38) as 

d^k ^ (45) 

A;2 + m2(r) 

Suppose rn?{T) > in (45). Then, in the Gaussian approximation, 

p{(j)^ = 0) = 7Ve-5^'/«/>.'^.>. (46) 

— * 

The momentum integration in (45) is cut off at = 0{T). It follows from (45) that 
^q{x) is correlated over a distance ^ = m(T)~^. 

The variance (44) can be written in terms of the Fourier transforms i{k) of the 
window function (normalised to /(O) — v) and 

G(k) ~ ^ , (47) 

k^ + m^{T) 



as 

< >=- l<fk \i{k)\^G{k). (4J 



Since the effect of / is to cut off the ffuctuations with wavelengths shorter than i?, 
where v = 0{R^), we can replace (48) by 

< 0^0^ >= 4 / Gik). (49) 

V J\k\<R-^ 

For a correlation volume v = 0{C,^), C, = m~^{T), it can be shown |Q that 

< 0„0„ >^ ATm{T), (50) 

where A ~ 10"^. Thus field ffuctuations on the scale of correlation volumes have rms 
value 

A0~(rm(T))i (51) 

[The prefector A is important in demonstrating that thermal ffuctuations are, in general, 
not likely to supplant bubble nucleation in first-order transitions [|19] but we shall not 
pursue this here.] 

3.3. The Ginzburg Regime. 

Now suppose that rn?{T) < 0. Still for a single field, 5*3 [$] can be written (cf.(38)) as 

SM = I d'x[^{\/^r + Ig'i^' - v\T)f]. (52) 

Consider the ffuctuations of the 'Higgs' field in one of minima, $ = 77 say, averaged 
over a correlation volume v = O {mJj^{T)) = 0{g^^r]{Ty^) for the Higgs field. As 
before, A<l> ~ {Tgri{T))2). On going close to the critical temperature Tc = 21] a point 
is reached at which A<l> ~ ri{T) . At this stage correlation volumes of the field can 
ffuctuate from the true vacuum at = ri(T) (or = —rj{T)) to the false vacuum 
(0 = 0) with significant probability. The temperature at which this happens is the 
Ginzburg temperature Tq 10 which, from above, is equally defined by 



A0(Tg) - viTc) (53) 

or 

gTc q Tg , . 

93 = = ~ 1 (54) 

r]{TG) niHyTG) 

or 

(1 - = 0{g\ (55) 

Equation (53) tells us how close we can get to the critical temperature before ffuctua- 
tions get too violent. The equivalent equation (54) says that the effective dimensionless 
coupling constant of the three-dimensional theory (which, from (36) and (38), is (73) 
becomes strong. Analysis of higher-loop diagrams would also show that it is at T = Tq 
that the one-loop approximation, with its mean-field behaviour, breaks down. Some 
care is needed on trying to define nearer to the critical point. As a three dimensional 
coupling constant it runs to an ultraviolet fixed point PD[ i.e. g vanishes with ri{T) (or 



miT)) at T = Tc. However, for our purposes we do not need to get closer to Tc than 
Tg and hope that the one- loop result remains qualitatively correct. 



In the symmetric phase the extension to a global 0{N) theory is straightforward. 
$ now becomes an N-component column vector. In the Gaussian approximation, the 
field distribution probabihties pt^ [$] become, from (35) 




(56) 



where Gab is diagonal, as 



Gab{x- X') = SabG{x- X'), 



(57) 



with G{x — x') given by (45). Root mean square fluctuations in |$| over a correlation- 
volume are 0{{Tm{T))2 , as before. 

In the broken symmetry phase the situation is more complicated because of the 
masslcssness of the Goldstone modes, In a Gaussian approximation Gab is no longer 
proportional to 5ab- However, the 0{N) radial field (Higgs) fluctuations must be such 
as to take |0| in a correlation- volume from a true vacuum at = rj{T) to the false 
vacuum at = with significant probability at a Ginzburg temperature Tq satisfying 
(55). For local gauge theories the scalar field is gauge- variant and we cannot define 
probabilities for arbitrary field configurations. Rather, we shall be interested in the 
probabilities p((5[0] = (l) for gauge- invariant Q. These will be developed as needed. 

4. STRINGS AND MONOPOLES NEAR 

So far our discussion of phase transitions has made no reference to the possible defects 
in the 0{N) theories. However a symmetry-breaking transition is implemented, the 
0{N) vector field does not switch uniformly in space (and time) from its symmetric 
vacuum value = to some value = 0c, \4>c\ — V on the vacuum manifold Ai = S^~^ 
which characterises the vacuum at large times. Rather, initially we expect to relax 
towards different elements 0c of Ai in different regions of space. What regions and how 
is the problem that we shall attempt to solve. Whatever, at some fairly early stage 
after the transition we expect a domain structure to have formed, such that within a 
single domain the 0{N) vector field is strongly correlated. It is between the domains 
that defects will be found. 

4.1. Global Strings. 

If we take a closed path through several domains the field takes a closed path on Ai. 
If Ai is not simply-conncctcd this path in Ai may be non-contractable. In that case 
the path in space has trapped a string (vortex) or strings. However, a non-contractable 
path is a 'large' path in field space. For there to be a reasonable chance for such a 
path to be executed the differences in the fields between adjacent domains must be 
substantial. That is, field fiuctuations must be large. The idea that defects are born 
out of large fluctuations of domains was flrst articulated by Kibble. 

We have just seen that one circumstance in which there are large fluctuations is 
close to a second-order phase transition, as characterised by the Ginzburg temperature. 
For global theories we know that we have the second-order transition for this to apply. 
We stress that this is not the only situation in which large fluctuations arise, and we 
shall turn to others later. If we knew more of the nature of phase transitions in the early 
universe we could determine whether Tq has any relevance. We don't know, but because 



it is the simplest option we examine it first, following an argument due to Kibble. [There 
are several mechanisms proposed by Kibble for the production of defects, according to 
the order of the transition, the timescales, etc.. See these proceedings. This is the 
simplest.] 

Specifically, we have seen that the Higgs field is correlated over a volume of 
characteristic size where ^ = m']^{T). From our earlier discussion, the first time 
that one can define such patches as the temperature falls below is at the Ginzburg 
temperature. In each patch the Higgs field is in the vacuum manifold but this in itself is 
not enough to determine the frequency of strings, since it is the field phases, rather than 
field magnitudes, that enforce their presence. The next step is to assume that the field 
phases are correlated on the same length scale ^ (even if they are Goldstone modes). 
In the Ginzburg regime the fluctuations are expected, from our previous analysis, to be 
so large that each patch will have a phase capable of differing by a large amount 0(1) 
from that of its neighbours. Thus, if it is further assumed that the phases in different 
correlation volumes are randomly distributed, we get large random jumps in phase as 
we cross domain boundaries. This is all that we need to show the presence of strings in 
large quantities. 

Instead of having to make these assumptions, they should be derivable for equilib- 
rium field theory at the Ginzburg temperature. We shall spend the rest of this chapter 
seeing to what extent this mechanism, the Kibble mechanism, is justified in this simple 
form, both for strings and for monopoles. 

By definition, if M. is not simply-connected then its first (fundamental) homotopy 
group ni(A^) is non-trivial. Since n„(S''") is trivial unless n = m, when n„(S'") = Z, 
in D = 3 dimensions there will be 0{N) strings only if = 2 i.e. 0(2) or U{1). The 
elements n E Z are the winding numbers of the strings, a change in the phase of the 
complex f/(l) field 

= 01 + i(f)2 (58) 

of 27m along the loop corresponding to a winding number n. For the field phase (f)/\4>\ 
to change by 27r (or a non-zero multiple of 2tt) requires that it cannot be well-defined 
everywhere since is continuous and single- valued. As a result must vanish somewhere 
in the loop. 0(2) or U{1) strings are thus characterised by lines of zeros of the field 
doublet, i.e. they are tubes of false vacuum, for which ~ 0, embedded in a true 
vacuum in which |0| ~ t]. 



The properties of global 0(2) strings are well-documented In the present 

context of thermal equilibrium they are the local axisymmetric instanton solutions to 
the field equations 

^ = -V^0. + i/0.(0^-r/^) = O. (59) 

The 0a in (59) have only spatial arguments 0^ = 0a(a;), since dimensional compactifica- 
tion has already been effected in our use of S3. More generally they are the solutions to 
the four-dimensional Euler-Lagrange equations for the action (11) i.e. solitons. How- 
ever, as written in (59), we have traded dynamical degrees of freedom for temperature- 
dependence of the parameters of the theory in a straightforward way. 

In cylindrical coordinates (r, 9, (p) a single string with winding number n, as solu- 
tion to (59), takes the form 



= r/(T)e^"7(r^7r/(T)), 



(60) 



where 

/(r) ^ CnT"", r ^ 0; /(r) ^ 1 - 0{r-'^) , r ^ oo. (61) 

As a result, we can deduce that 

a) the string thickness a(r) is 0(C(T)) = 0{m^(^)) = 0{g-\(^)-^), the Higgs 
field correlation length at temperature T. 

b) strings with winding number n > 1 (which can be thought of as lines of multiple 
zeros) are unstable, preferentially splitting into strings for which n = 1. 

c) the energy per unit length diverges logarithmically, cut off by the presence 
of nearest neighbour strings in any network. Logarithms apart, the energy/length is 
0{r]'^{T)), vanishing as we appproach the transition. 

d) the net winding number of strings through a surface S is given as the line 
integral 

Ns = ^i dl.da, (62) 
2ti J as 

where we have adopted the radial/ angular field decomposition 

= pe^". (63) 
In the complex field notation of (58), Ns of (62) can be reexpressed as 

More usefully, as a surface integral, it is equivalent to 
or, in terms of p and a 

TT JS p^ 

To see that fiuctuations near the phase transition are indeed capable of creating 
vortices we need to calculate the probability that the winding number of the field 
through a loop dS bounding a surface ^S" be n. In practice it is more convenient to 
evaluate the related quantity 



No = —i dl. 



27r Ids 'v{T)\(f)\ 

I d~s'.(dpAda) (67) 

s 



2tt7]{T) Js 

For a large loop dS the difference between Ns and Ns (not integer) is vanishingly small 
if no vortices pass close to dS, and Ns remains a good indicator of vortex production. 
In the Gaussian approximation the probability that Ns{t) takes the value n is 

p(Ns = n) = exp{-in7 < NgNg >}. (68) 

On decomposing the radial mode as p = r]{T) + h for Higgs field h and defining the 
Goldstone mode g hy g — ri{T)a, from (67) it follows that 

< NsNs >= {^^^^y J J^dS'dS"<{dh'Adg'){dh"Adg")>. (69) 



The primes (doubleprimes) denote fields in the infinitesimal areas dS', dS" of S respec- 
tively. For economy of notation we have not made ths scalar products explicit. Without 
loss of generality we take S in the 1-2 plane, whence 

< NsNs >= i^^^^fj j,.^,/^'^^" < ^^h'^ih"^Jg'^Jg" - dih'd^h"d^g'd,g" > . 

^70) 

It is convenient to refine our notation further, decomposing space-time as a; = {t, x) ~ 
{t, xl, xt) where xl = {xi, X2) denotes the co-ordinates of and xt — x^ the transverse 
direction to S. Similarly, we separate 4-momentum p as p = {E,pl,pt)- ■ 

Let Ghix'-x") =< h{x')h{x") >, Gg{x'-x") =< g{x')g{x") > be the Higgs field 
and Goldstonc mode correlation functions respectively as read off from 6*3. As a first 
step we ignore correlations between Higgs and Goldstone fields. That is, we retain only 
the disconnected parts of < Ns^s >■ Eqn. (70) then simplifies to 

< NsNs >= (tt^T^)' / / "^^'dS" [< dih'dih" X djg'djg" > 
27rf]^{I ) J Ji,j=i.2 

- < dih'djh" >< djg'dig" >] (71) 

which can be written as 

(72) 

In (72) I{pl) is the Fourier transform of the window function I{xl) of the surface 5" 
(i.e. I{xl) — 1 ii Xl & S, otherwise zero). 

We coarse-grain in the transverse and longitudinal directions by imposing a cut-off 
in three-momenta at < A = for some /, as before. Thus Ns is now understood 
as the average value over a closed set of correlation-volume 'beads' through which dS 
runs hke a necklace. 

For large loops dS, I{qL) —^Ql), enabling us to write 

< NsNs >= (73) 

{^rJ(tqTja'pG,{p + qT,t)G,{p,t) j<fqL\i{qLm{pLf{qLf-{pL.qLf] 

By p + we mean {pl,Pt + qr)- The dependence on the contour dS is contained in 
the final integral 

J ^ jd-'qLmLmipLfiqLf-ipL.qL?] (74) 

~ TTpi fdqLql\I{qL)\' (75) 

If this is evaluated for a circular loop of radius R, we find 

J^plO{27rR/l), (76) 

as we might have anticipated. The rms winding number behaves with path length as 
An oc J'2 = 0(L2 ) , where L is the number of steps of length /. This is consistent with 
the coarse-grained field phases of different volumes v being randomly distributed. 

The final step is to relate the magnitude of the fluctuations in winding number 
Ns to the magnitude of the radial (Higgs) held fluctuations and angular (Goldstone) 



field fluctuations. To estimate the fluctuations we a) neglect the positivity of p and the 
Jacobian from the non-linear transformation (43) and b) the non-singlevaluedness of a. 
While valid for small fluctuations about the global minima this can only be approx- 
imate for large fluctuations. With these provisos, at temperature T the equilibrium 
propagators are time- independent, dominantly the free- field propagators 

+ mjj[T) 

GM = ^ (78) 

If we take p| = 2p^/3 in the integral then, up to numerical factors, we can approximate 
(4.16) as 

< NsNs >^ J (t;^) mH{T)T f dpGh{p), (79) 

VzTTCr^/ J\p\<mH{T) 

where we have coarse-grained to the Higgs correlation length ^ = mjl{T). [We have 
further assumed that the qt integration can be approximated by setting qx to zero in 
the integrand. Qualitatively this is a reasonable simplification]. The integral in (79) is 
essentially the integral (49). The end result is that after substitution, 

<iV.iVs>=o((^^)yO(L) (80) 

where L = RmniT) is the length of the path in units of ^. Equivalently, on using our 
previous results for equilibrium 

<NsNs> = Qi I <^A-> y (81) 




O f7^ \0{L) (82) 



Thus, as anticipated by Kibble, the phase fluctuations are indeed scaled by the Higgs 
fluctuations. Further, in the Ginzburg regime, when correlation-volumes of the Higgs 
field can fluctuate to the false vacuum with significant probability, the fluctuations in 
field phase on the same distance scale are of order unity and, from the 0{L) term, are 
distributed randomly. 

Even though fluctuations are large enough to produce non-trivial winding number 
there may still be some doubt that this is directly related to strings. However, strings 
intersect sufficiently infrequently, at low density at least, for energy to be approximately 
proportionately to length, with a calculable entropy density 

Sst = 0{g7^{T)) = 0{mH{T)). (83) 

The effective string tension, whose vanishing signals string proliferation at a Shock- 



ley/Hagedorn phase transition, is [ p^ 

(Jeff = cr -Tsst, (84) 



where a = 0{7f{T)) is the string energy/unit length. On inspection it is seen that (Te// 
also vanishes at the same Ginzburg temperature Tq of (55) for which the fluctuations 



in Ns become large. Since strings produced from random field phases approximately 
execute random walks, the majority of string lies in 'infinite' strings, rather than small 
loops of length O (^) . [The probability of a self-avoiding walk not forming a loop is ap- 
proximately seventy percent]. Thus, at Tq it costs nothing to produce such macroscopic 
defects from local fluctuations. 



4.2. Global Monopoles. 

If, instead of taking a closed path through the domains after the transition, we take 
a closed surface through them, the 0{N) field vector on the closed surface executes a 
closed surface in JH = S^~^. If, in turn, this closed surface is non-contractable, then 
point defects (monopoles) can be trapped within it. As with strings, the formation 
of monopoles requires large field fluctuations. The noncontractability of surfaces in Ai 
requires the second homotopy group 112(5'^"^) be non-trivial. In D = 3 dimensions this 
only happens when N = 3. Thus, at temperature T, a global 0(3) theory possesses 
global monopoles (balls of false vacuum), for which the simplest solution (winding 



number n = 1) is of the form ^1 



US)=r]{T)higr]{T)r)^, (85) 

r 

where r now denotes radial distance, and h is to be determined from the Euler-Lagrange 
equations. We are less interested in monopoles than strings for cosmological reasons. 
Sufficient to say that 

a) the monopole diameter a(T) = 0(^(T)) = 0{g^^r]^^{T)) = 0{m]j^), the Higgs 
correlation length 

b) the monopole energy within a sphere of radius r diverges linearly with r. 
Nonetheless, as with strings, the presence of nearby defects cuts off the energy di- 
vergence. Of greater interest, 

c) monopoles with winding number n > 1 are unstable 

d) the net winding number of the monopoles in a volume v is 

In calculating the variance in winding number in a volume v a similar analysis can 
be performed to that for strings. Instead of Ny, we calculate the variance of its relative 



Ny = — ^ / d'x '^"77°, ^^ (87) 
87r J r]'\T) ^ ' 

Then, in an approximation that ignores the connected correlation functions i.e. re- 
places the six-point correlation function by the product of three independent two-point 
functions, we find 



< NyNy >= 0{{RmH{T)Y){g^T/m{T)f, (88) 

where v = 0{R^). The behaviour, counting the number of correlation patches on 
the surface of v is generic, interpreted as refiecting the random nature of the 'angu- 
lar' distribution of the Goldstone modes. The remaining factor, derived here in an 
approximation in which free- field propagators have been assumed dominant, links the 



correlations of the Higgs field to those of the Goldstone modes. As with strings, in the 
Ginzburg regime, fluctuations are sufficiently strong for field phase differences between 
Higgs correlation volumes to be random and of order unity, the prerequisite for a large 
monopole presence. 

So far, everything accords with the picture proposed by Kibble in that the fluc- 
tuations at the Ginzburg temperature are large enough to contain defects with the 
behaviour that we would have expected. However, the approximations that we have 
been obliged to make are more severe near the transition, where we most need reli- 
able results. In particular, the neglect of connected correlation functions means that 
the variance < NN > does not know whether defects are present. We conclude these 
observations on fluctuations by examining the strings of the local U{1) gauge theory. 
Here the picture is complicated by gauge properties, which have lead some authors 
to doubt the Kibble mechanism as presented so far. However, as will be seen, in other 
ways they permit more justifiable calculations. 



4.3. Local Strings 

The arguments about symmetry breaking and the existence of defects relied only on the 
group structure of the vacuum manifold. Gauging the global symmetry does not effect 
these arguments as the gauge fields are zero in the vacuum. The precise dynamics of 
string formation and interactions will of course be rather different in the local case. For 
instance, it might be that in gauge theories the density of strings formed at a transition 
is far too low to be of much interest. 

One obvious hurdle which comes with the introduction of gauge symmetry is the 
fact that parameters which were physical in the global case become gauge dependent 
and hence unphysical in the case of local symmetry. In our example of strings in the 
0(2) theory, the different vacua are labeled by different values of the phase a of the 
field (when thought of as a single complex f/(l) field). However, a is gauge dependent. 
Indeed in the unitary gauge it is removed altogether! It is the degree of freedom which 
is eaten by the gauge fields as they become massive. The question is, does the move 
from global to local symmetry effect the basic Kibble mechanism arguments? 

One picture that might be used is to fix the gauge to be something fairly straight- 
forward (the so called unitary gauges are well known to be difficult and so we exclude 
them). With the gauge dependence dealt with provisionally, we can use unphysical fields 
as convenient intermediate tools provided we end up studying a physical quantity. 

So the arguments now proceed as they did in the global case. Ns of (4.5) is 
equally valid for the local U{1) gauge theory. One can then consider following the 
phase as one moves along a closed loop dS in space, keeping track of the number of 
times the vacuum manifold is traversed. The phase must return to the same point 
on the manifold so that / a = 2ttNs where Ns is the winding number, an integer. 
As the winding number is a physical quantity it does make sense to discuss the line 
integral of the phase even though the phase is not itself physical. In moving round the 
loop of length R, one passes through R/^h correlated patches in the Higgs field, for 
which = m^. Provided the fixed, but gauge- dependent, correlation length for the 
field phase a is, on average, scaled by and the phases are randomly distributed in 
the domains, then from a random walk analysis one expects that (NsNs) = R/C,h, as 
before. 



One might remain uneasy about relying so heavily on an unphysical parameter, 
and there have been attempts to show that this picture is wrong. We therefore turn 
to a second approach which can be used in the 0(2) or U(l) case. For this we to look 
at the classical solutions which minimise the compactified action 5*3 in which all heavy 
modes have been integrated over and which correspond to strings. These are rather 
more complicated than their global counterparts. Nonetheless, there are well known 



vortex solutions to these equations |^ in which there is an infinite straight (necessarily) 
static string with winding number n. 

The key observation is that if the scalar field is to be in the vacuum at distances 
large compared to the string size (set by correlation lengths) then 

eA^' = d^'a . (89) 

Indeed on studying the exact solutions we find that in cylindrical coordinates (r, 9, z) 
with the string along the z axis has as asymptotic solutions at large distances 

Ae = -- cy \ e-^(^)'- , a = ne, |0| = r,{T) - cu . .-^(^^ 

(90) 

with other components of A^^ being zero. The effective Higgs mass mniT) is gr]{T) as 
before, and the effective gauge field mass is my = er](T). This satisfies (^) up to small 
corrections. The idea then is to use the gauge field to track the Higgs phase and hence 
count the strings passing through some loop. There is a critical value of e/g below 
which strings with n > 1 are unstable and we assume this to be the case, since it is 
approximately correlated to the existence of a second-order transition. 

To calculate the fiuctations in the numbers of strings Ns which pass through a 
surface S we can now use (p9|). The variance in Ns is then expressed completely in 
terms of a two-point correlation of the gauge fields, 

2 



dxi f dx'j.{Ai{t,x)Aj{t,x')), (91) 



47r2 

4:71^ Ids ' fdS 

rather than the four-point scalar field correlation function (4.12) of the global theory. 
Thus we just need to know the equal-time correlation function of the A fields. There are 
only two symmetric tensors in the spatial indices which can be made in the covariant 
and gauges, namely 5*-' , k'^k^ . The Ward or BRS identities and the Lorentz structure 
at finite temperature then tell us that the most general form for the (time-dependent) 
correlation functions in covariant and i?^ gauges is 

tG'^{t-t',x-x') = {TA\t,x)A^{t',x')) (92) 

+ 00 „ 

^ d^ke^'^^''^'-''^-^-^^'^'^\G'^{k)} (93) 



n=—oo 



1 P k^ 



where the energy ko = 27rm//3. The / function varies with the gauge chosen while the 
other and 11^ terms correspond to the physical modes of the photon, two transverse 
(magnetic) and a single longitudinal (electric) mode. 

In the case of global strings it was sufficient to take the n = equal-time term 
in the sum, which could be identified with the compactified action S3, in terms of 
which fiuctuation probabilities were given. We can do this here, even though is 
nonphysical, but will not. The reason is that, later, we need to know some aspects of 
its time behaviour. 

From (pri), (P^, ( P^ we see that we require two sorts of integration 

c/xiC-^^-^ / dxiC-'^-'^h. (95) 

dS JdS 

Since the latter is zero we only pick up the term in the propagator with the This is 
not surprising as we are really looking at the magnetic field correlations. The transverse 
correlation function with 5ij is associated with magnetic fields, its self-energy TIt giving 
information about magnetic screening. Thus we have that 



e2 
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dxi e 

dS 



(96) 



Note that this form is independent of the gauge chosen. 

In order to restore a temperature dependence to 77 we calculate Il.T{k) by resum- 
ming an infinite set of diagrams in the consistent manner of Braaten and Pisarski. The 
major effect will be that IIt provides the required shift in the expectation value 77. 
However there are further corrections due to the presence of strings in the background. 
This means that 11^ is a non-trivial function of momentum leading to both cuts and 
poles in the propagator in the complex energy plane. However, for a low density of 
strings, we know that the solution is almost the same as being in the vacuum state so 
as a first guess we need not consider corrections to Hy due to the strings' presence. 
Thus we work with a photon mass of the form 

eW + ^T{k) = e'viTr + <ag = M^ (97) 

where we have explicitly included the self-energy effects which give a temperature de- 
pendent 7] and left other effects in the m"^^^ term. The mass rrimag is the magnetic 
screening mass that, in general, would be expected to be present. For the special case 
of a U{1) theory rrimag = if the string background is ignored. We can now do the 
energy sum to give 



{NsNs) = ^Jd'k 



32tt^ J 00 exp{(3uj} 



dxi e 

dS 



+ (T=0 terms) (98) 



where cu^ = fc^ + M^. 

If we now limit ourselves to counting N within a circular loop dS of radius R, we 
find that that in terms of the Bessel function 



/ dxi i dx'i e-^^-(^-^') = i?2 dedd' cos{e - 9') exp{-iRk.{cos{e) - cos{e'))} 

JdS JdS J-n 

= A-n'^R\h{kR)\^ (99) 
r dd exp{-z2. cos(^) + inO} (100) 

J — TX 



g — i7i7r/2 



27r 



and so we have that 



p'i T?2 ^oo POO _ _ 1 1 _ 

(NsNs) = ^/ dkj d~k27rk , \Ji(kR)\'' (101) 

— * — — * 

where kz is the component of k perpendicular to the loop and k the size of the k in the 
plane of the loop, so that k^ = k"^ + k^. 

We now look at very large loops, and high temperatures T ^ m ^ R~^. We can 
then use the asymptotic expansion for Ji 

(M^))'^ ~ —[1 + cos(2z - 37r/2).] (102) 

TTZ 

Thus we find that the largest contribution comes from a region where m < k < T so 
that 

(NsNs) ^ ^^Tdu; \ (103) 

which using the high temperature limit gives 

(TV.TV.) ==—!„(-) (104) 

As we observed, this result could have been obtained directly from the n = term 
alone in the mode sum, which corresponds to replacing exp{/3a;} — 1 by (5uj in the 
denominators, and remembering that compactification provides a cutoff at the thermal 
wavelength \k\=0{T). 

Superficially < NsNs > of (4.47) looks very different from its counterpart (4.23), 
apart form its linear dependence on the loop perimeter R. However, if we choose 
^h{T) = mff^{T) as the natural length scale, then (4.47) can be written 

(105) 

where we have set rrimag to zero. Assuming that e/g ^ 1, compatible with a second- 
order transition, then (up to logarithms), the scale of < NsNs > is set by = g'^T /itlh- 
In the Ginzburg regime, where g^ ~ 1, the variance in Ns is 0{R/^h), with coefficient 
0{ln{l/eg)). Qualitatively, this is as before, provided eg is not extremely small. That 
is, once in the Ginzburg regime, fluctuations are large enough to contain strings with 
high density. 

However one question that remains unanswered is exactly how does (NsNs) know 
if there are strings around or not? The result is only sensitive to the presence of strings 
through the \4>\ dependent mass term and we know that, for dilute strings at least, the 
correction to due to strings being present is small. This can be made more rigorous 
by carefully calculating the photon propagator in a nontrivial background which is a 
minimum of the classical action. If one leaves the terms cubic or higher in the action 
as perturbations to be treated in the usual way, it is easy to see, in the first instance, 
that the photon propagator has the same form as before but where the |0| in the mass 
term is replaced by the classical scalar field solution. In the case of dilute strings this 
whole procedure is justified as the deviations this produces from calculations in a the 
usual fiat background of the true vacuum are exponentially suppressed. The end result 




is that, although the fluctuations in the field are reduced, it seems that the result (104) 
remains qualitativly correct. The concerns that a local theory would possess strings 
with a lower density than a global theory is unfounded in this particular instance. 

We should stress that the 0{N) theory, with its Z-strings, may be special in some 
regards. We have little idea as to the nature of the symmetry-breaking at GUT scales, 
but it is equally likely that the relevant homotopy group is Z2, rather than Z. Since 
the net winding number of Z2 strings is 0, 1 the previous analysis is inappropriate and 
we don't know where to begin. 

Finally, the local O (3) theory permits t'Hooft-Polyakov |^ monopoles. Unfortu- 
nately, their winding number variance cannot be expressed simply in terms of gauge 
fields and we shall not consider them. Since, in the early universe, monopoles are an 
embarrassment, this is no great loss. 



5. TIME SCALES 

In most of the calculations discussed hitherto, it is assumed that the system is in thermal 
equilibrium, enabling us to use the imaginary-time formalism. At best, in early universe 
calculations, equilibrium can only be achieved in some intermediate period. Prior to 
this, as we noted in the introduction, the universe is expanding too rapidly for the 
particles to interact sufficiently to equilibriate. Subsequent to this particle species, or 
defects, will freeze out for chemical or thermal reasons. It is worthwhile spending a 
few moments to see whether this temporary lull can occur for a temperature range 
encompassing the critical temperature appropriate for an 0{N) GUT transition (or 
that of any other symmetry group). If this is not possible our work of the previous 
sections is irrelevant for the very early universe. 

Let us assume equilibrium, for which it is relatively easy to estimate the time 
scales. We can then check for consistency. The time scale associated with the rate 
at which the universe is expanding Texp is simple to quantify. It is the inverse Hubble 
constant 

where we have used the Friedmann equation to link the Hubble constant to Newton's 
constant. Here (7* is the effective number of degrees of freedom |^ (total number 
of relativistic degrees of freedom with a few modifying factors). It is around 100 for 
T > IGeV in the standard model. The time scales for particle interactions are also easy 
to calculate. Looking at the imaginary part of the self-energy tells us the time taken for 
a particle of a given frequency to interact with any possible other real-physical particles. 
For particles stationary in the heat-bath we find at high temperatures (T larger than 
all other particle physics scales) 

The first expression is for a gauge boson in a pure SU(N) theory. To obtain the precise 
numerical factor requires considerable effort but this is essential if one wants to have a 
gauge invariant and therefore believable answer [ES]. In the early universe contributions 



from fermions and other matter will decrease this expression further. The decay rate 
for a fermion in a gauge theory takes a similar form except that the N is replaced by 



group theory factors depending on the representation of the fermions involved ||29|. One 
should note that there are unresolved problems with the thermalisation rates in gauge 
theories of particles moving with respect to the heat bath because of infrared problems 
associated with the behaviour of magnetic fields in a plasma . 



The latter time scale in ( |107|) is for self-interacting scalar s |31] with interaction 
(7^0^. The reason for using this notation for the scalar coupling is that for most quan- 
tities, it is when e ~ (7 that the two types of interactions give similar-sized effects. For 
instance the thermal mass corrections are 5m ~ eT, gT . One of the exceptions to this 
simple rule are the thermalisation rates quoted here. This is also marks the boundary 
between first and second order transitions as mentioned at the end of section 2. 

These time scales are an estimate of the time in which a particle traveling through 
a plasma starts to lose its energy to the heat bath and so becomes absorbed as just 
another part of the thermal background ||3^. They only tell us about the start of the 



energy loss process as the result is to be interpreted through linear response theory in 
which one considers only small or short time deviations from equilibrium. Essentially 
these results are telling us that the calculations which produced them become invalid 
and inconsistent at this time scale. To follow the process of coming into equilibrium 
properly requires that one goes beyond the usual linear response analysis, which is very 
difficult to do. 

A simpler approach is to consider a single interaction. The crossection goes as 
a ~ e^/s,g^/s where s is the characteristic centre of mass energy squared scale which 
here will be T^. The many-body nature of the problem can then be taken account. 
First the density of particles per unit volume in equilibrium is n ~ and second their 
speed, is 0(1). An estimate of interaction time is then 

111 , , 

Note that this simple argument does not give the correct behaviour for the scalar case 
as given in (|107|) . 

For equilibrium we require that 

rexp>r ^T</Mpi, e^Mpi, (109) 

where Mpi is the Plank mass. Although qualitatively correct, the prefactors are impor- 
tant. Assuming the Standard Model degrees of freedom, the inequalities become 

T « e'^Pi N ^ ^ iQifiGeV, (110) 

1.6gl' 11-1 

T < 2.5g^ X IQi^GeV (111) 

Thus gauge particles are in equilibrium at GUT scales and below, at least until the high 
temperature forms for the particle interactions (|107|) and ( |108|) cease to be valid. The 
second inequality ( |111| ) is the result for a pure 0(1) self-interacting scalar, suggesting 
that they are not usually in equilibrium at GUT scales. 

In fact, these temperature upper bounds may be an underestimates. There has 
been a suggestion that the characteristic relaxation time of the plasma is actually ||3^ 



l/(e^T) rather than ( |107| ). We are unsure about this. Apart from using more sophis- 
ticted formulae for particle interaction times, valid over a larger range of temperatures, 
the simple analysis given above does not take phase transitions into account. Second 



order processes show critical slowing down. Rather, the purpose of the calculation is to 
see whether thermal boundary conditions are appropriate. 

The above calculation shows that most species are in equilibrium below GUT 
temperatures. It is however possible for a species to fall out of equilibrium, to 'freeze' 
out at a later stage. This process is extremely relevant to defects. However the principle 
is most easily illustrated with the photons of the microwave background. 

At recombination, the total crossection for all photon interactions falls dramat- 
ically over a short time period. This essentially leaves a background of thermal non- 
interacting photons of temperature T ~ 0(eV). The number and energy distributions 
of the photons are not changed by particle interactions to a very good approximation. 
The photons still interact gravitationally but this simply involves a dilution of the pho- 
ton density due to the usual stretching of the universe. It is easy to check that the 
photons keep a black-body spectrum and that the expansion of the universe causes 
their effective temperature to fall from around lO'^K to the present 2.7K. 

More generally, the freeze-out can be seen by using formulae for particle interaction 
time scales appropriate at lower temperatures. Essentially, the time scale goes from 
e~^T~^ to e~^A"T~'^~^ where A is some other particle physics scale relevant to the 
particular species. If n > 1 then we get a lower bound on the temperatures until which 
a given species remains in equilibrium 



This is the freeze-out temperature of the particle species. 

This is a particularly simple example because the photons become essentially non- 
interacting (except for gravity). There are however different shades of equilibrium and 
corresponding freeze out as species fall out of 'equilibrium'. What we have in mind here 
is the difference between thermal and chemical equilibrium. So far we have concentrated 
on thermal equilibrium. By this we mean the situation where the energy distribution 
of a species is being maintained in a Bose-Einstein or Fermi-Dirac distribution. This 
sort of equilibrium will be maintained by interactions with any other species already in 
equilibrium and this includes self-interactions. So to test for thermal equilibrium one 
need only check that the time scale associated with the fastest interaction process is 
shorter than the expansion rate of the universe. 

Chemical equilibrium is concerned with the maintenance of equilibrium particle 
numbers. Each individual interactions mixes only certain types of particle species to- 
gether. The quark-gluon plasma provides an excellent example which has been studied 
extensively becuase of the application to terrestial relativistic heavy-ion collision exper- 
iments |3^. If one heats up nucleons to obtain a quark-gluon plasma at a few hundred 
MeV, as is done in relativistic heavy ion collisions, one starts with an initial domination 
of up and down quarks plus gluons. The gluons thermalise about three times quicker 
than the quarks so initially the plasma is pure glue. While not strictly in equilibrium, 
the gluons maintain a Bose-Einstein distribution. The up and down quarks then come 
into equilibrium in which there are significant numbers of both up and down quarks 
and anti-quarks. There is an imbalance between the quark and anti-quark numbers 
reflecting the initial baryon number. Strange quark- ant i- quark pairs can be created 
as easily as up and down quark pairs as at these temperatures the mass differences 
are being negligible. It will, however, take a considerable length of time for the im- 
balance between quarks and anti-quarks to to be felt in the strange quark densities as 



T 




(112) 



the stangeness changing processes are very slow. The time scale of such strangeness 
changing interactions, is much longer than other quark interactions Tq. In the rela- 
tivistic heavy-ion collision experiments this leads to a situation where on time scales t 
where < t < the plasma can have large patches where the energy density of the 
total number of quarks follows a Fermi-Dirac distribution. The net baryon number is 
carried almost entirely by the up and down quarks whereas in the long time limit it 
should be spread evenly between all the quark flavours which are essentially massless at 
these temperatures. This can be modeled by having the number of up and down quarks 
approximately conserved and the number of strange quarks also conserved. This means 
that one would have different chemical potentials associated with the total number of 
up plus down quarks and for strange quarks instead of one associated with the total 
number of all types of quark. These chemical potentials would slowly change on time 
scales of order r^. 

In the present context of defects, the interesting question is not if the kinetic energy 
in the defects follows a thermal distribution but what is the density of defects plus anti- 
defects. This is what is important for the gravitational effects, for instance. Thus we 
are asking about the chemical equilibrium of defects, not the thermal equilibrium. The 
question then is what is the time scale for processes which can change the net amount 
of string passing through a loop or the net number of monopoles in a given volume. 

In the usual picture the universe settles down into equilibrium after a phase tran- 
sition and one then finds that the defects are frozen in. That is to say, it is assumed 
that fluctuations are no longer strong enough to change the numbers of zeros in the field 
which represent positions of monopoles or strings in the 0(3) and 0(2) cases considered 
here. To verify this we are therefore interested in trying to calculate the fluctuations in 
the numbers of defects, and more especially the time-scale over which such fluctuations 
occur. This is what we have to compare to the expansion rate of the universe. 

At first sight it might seem as if it is very unlikely that a fluctuation could change 
the total string or monopole number. They are, after all, topological charges and to 
change the total winding number of the phase at infinity for instance would require large 
changes in the field at all points in space. However, what is relevant is the total defect 
density and not the defect minus anti-defect density, which is the conserved number. 
Thus we would like to try and count how often fluctuations create loops of string or 
monopole anti-monopole pairs for instance. 

This is not an easy task. One simple approach would be to take the estimate of 
the fluctuations in the numbers of local strings {NN) and now ask how does this change 
in time i.e. to calculate 



where N{t) is the number of strings minus anti- strings passing through the same 
loop in space but at different times t. We have assumed thermal equilibrium so that 
{N{t)N{t)) = (A^(O)A^(O)). It is easy to upgrade the calculation of {N{0)N{0)) made 



{[N{t) - iV(0)]2) = 2(iV(0)iV(0)) - 2{N{t)Nm, 



(113) 



above. From ( pT]) we have 




(114) 



This can be expressed in terms of the propagator as before at the expense of one extra 
exponential in time 



(iV(t)iV(O)) 
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In doing the energy sum we can convert into a contour integral in the usual way. With 
the same approximations made on the propagator, the energy sum leaves contributions 
coming from the residues of the poles of the integrand which are at /cq = ±1^ where 
uj = {k + M^)^/^ as before. We merely pick up an extra factor of e**"^ in the remaining 
integrals. 

Putting all this together the generalisation of (jOT 



IS 



([iv(t)-iv(o)n 



d^k 



1 2sin2(cjt/2) 



327r^ J UJ exp{/?ti;} 
+ (T=0 terms) 



-ik.x 



(116) 



Without going into details it is clear that this becomes large when ~ cu, where uj is 
effectively limited to values between M ^ eT and T. Unfortunately, it is not clear in 
such a simple calculation if this if the time scale for fluctuations due to loops of string 
being created around the boundary of the loop (any inside the loop will contribute 
nothing as the loop measures string-anti-string numbers only). It could just be picking 
up the inevitable motion of the strings in and out of the loop. In short, it will take a 
more sophisticated approach that accomodates a string background, to try and get a 
measure of the time scale of defect number fluctuations. Until that time, the picture 
of defects frozen in shortly after the phase transition can not be confirmed by analytic 
calculations. 



6. NONEQUILIBRIUM BEHAVIOUR 

We have seen that equilibrium methods are not powerful enough to demonstrate how 
defects formed at the Ginzburg temperature can freeze in after the transition to survive 
as topologically stable (macroscopic) entities. For this we need nonequilibrium methods. 
As yet they too are poorly developed to deal with defect formation at the Ginzburg 
temperature. However, in one sense this may not be the issue. One reason why so much 
thought has gone into fluctuations at the Ginzburg temperature is that it was necessary 
to demonstrate at least one circumstance in which fluctuations could produce strings 
and monopoles in large numbers for defects to be viable, in principle. Now that we 
have done this we can move to more compelling pictures of the way in which the initial 
density of topological defects is flxed. In practice, while the domain mechanism outlined 
initially is almost certainly correct, it ccould well be that the Ginzberg temperature is 
relevant to nothing other than a thermally produced population of defects. 

We saw earler that an expanding universe is driving the system very fast. As 
the system is driven from some initial state there comes a point when the rate of 
change is too fast for the evolution of the fleld to keep up [0] [This is true even for 
very slowly driven transitions, but for appropriate couplings the out-of-equilibrium 
behaviour occurs more closely to the transition, in time and temperature, than the 



Ginzburg regime. See Kibble, these proceedings]. The transition may now be viewed as 
a quench and it is not clear that either temperature or free energy mean anything at all. 
At this point any defects within the field are assumed to be frozen in until the transition 
is complete. Upon completion, the field will try to establish thermal equilibrium. At 
sufficiently low temperatures, however, the return to equilibrium by thermal processes 
will be so slow that the evolution of the initial defect density thus produced will be 
almost entirely by interactions between the defects themselves. 

Thus, in this scenario, the vortices cease to be produced, not at the Ginzberg 
temperature, but when the scalar fields go out of equilibrium The relevant scale which 
determines the defect density is the coherence length, ^{t) at this time, and for some 
time onwards, rather than the coherence length at the Ginzberg temperature. The 
initial number of vortices produced during a phase transition would then be expected 
to be roughly the number when the scalar field first goes out of equilibrium. The next 
part of these notes will go some way in justifying this expectation. 

Gauge fields are too complicated to handle out of equilibrium. We shall return 
to the global 0{N) theory in D = 3 dimensions, with its monopoles for = 3 and its 
strings for N = 2. From the viewpoint above the transition is realised by the changing 
environment inducing an explicit time-dependence in the field parameters. Although 
we have the early universe in mind, we remain as simple as possible, in flat space-time 
with the 0-field action of before:- 

S[<l>] = I - lm'it)<l>l - Ig'itMr). (117) 

but for the t-dependence of m?{t) and g'^{t) , which is assumed given. 

We wish to calculate the evolution of the defect density during the fall from the 
false vacuum to the true vacuum after a rapid quench from an initial state. The simplest 
assumption, made here, is that the symmetry breaking occurs at time t = to, with the 
sign of m?{t) changing from positive to negative at to- Further, after some short period 
At = t — to > 0, m?{t) and A(t) have relaxed to their final values, denoted by and 
respectively. The field begins to respond to the symmetry-breaking at t = to but we 
assume that its response time is greater than At, again ignoring any mode dependence. 

To follow the evolution of the defect density during the fall off the hill towards 
the true vvacuum involves two problems. The first is how to follow the evolution of the 
quantum field, the second is how to count the defects. We take these in turn. 

6.1. The Closed Timepath Approach to Non-Equilibrium 

During a rapid transition the dynamics of the quantum field is intrinsically non- 
equilibrium. The normal techniques of equilibrium thermal field theory are therefore 
inapplicable. Out of equilibrium, one typically proceeds using a functional Schrodinger 
equation or using the closed time path formalism [^. Here, we employ the latter, 
following closely the work of Boyanovsky, de Vega and coauthors 0, |^ . 

Take t = to as our starting time. Suppose that, at this time, the system is in a 
pure state, in which the measurement of would give ^o{x). That is:- 



0(to,x)|<l>o,to) = $o|$o,to)- 



(118) 



The probability pt^l^f] that, at time tj > to, the measurement of (p will give the value 
$/ is Ptfi^f] = |c/oP, where:- 

c/o= / V(j)exp\iSm, (119) 

in which = Jl^i T^'Pa and spatial labels have been suppressed. It follows that pt_^. [$ f] 
can be written in the closed time-path form:- 

PtA^f] = / Vct>+Vct>. exp z 5[</)+] - . (120) 

Instead of separately integrating 0± along the time paths to < t < tf, the integral can 
be interpreted as time-ordering of a field along the closed path C+ © C_ as shown in 
figure where = 0+ on C+ and = 0_ on C_. It is convenient to extend the contour 
from tf to t = oc. Either 0+ or 0_ is an equally good candidate for the physical field, 
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Figure 3. Closed Time Path 

but we choose 0+:- With this choice and suitable normalisation, pt^ becomes:- 

Pt.l^f] = I P0+P0-5[0+(t) exp|z(5[0+] -5[0_]) j, (121) 

where 5[(t)^{t) — $j] is a delta functional, imposing the constraint 0_,.(t,x) = $/(a;) for 
each X. 

The choice of a pure state at time to is too simple to be of any use. The one 
fixed condition is that we begin in a symmetric state with (0) = at time t = to- 
Otherwise, our ignorance is parametrised in the probability distribution that at time 
to, (j)(to,x) = $(x). If we allow for an initial probability distribution -Pto[$] |0 then 
Ptj^l^f] is generalised to:- 

Pt,[<^f] = jvmM /^^^^^^^^P0+P0-5[0+(t)-$/]exp|z(5[0+]-5[0_])}. (122) 

At this stage, we have to begin to make approximations. To make any progress, we 
shall find that it is necessary that be Gaussian. For this to be so -Pto[$] must 

be Gaussian also, with zero mean. For all our caveats about the difficulty of defining 
temperature we don't know where to begin unless we assume that $ is Boltzmann 
distributed at time to at an effective temperature of Tq = (3q^ according to a quadratic 
Hamiltonian -ffo['^']- That is:- 

PtM = ($,to|e-'^^»|<I',to) = / P03exp(^^o[03]|, (123) 



for a corresponding action 5*0 [^s], in which 03 is taken to be periodic in imaginary time 
with period /Sq. We take 5*0 [^3] to be quadratic in the 0{N) vector 03 as:- 



(124) 



We stress that ttiq and /3o parametrise our uncertainty in the initial conditions. The 
choice Po ^ 00 corresponds to choosing the pt [$] to be determined by the ground state 
functional of Ho, for example. Whatever, the effect is to give an action 5*3 [0] in which 
we are in thermal equilibrium for t < to during which period the mass m{t) takes the 
constant value mo and, by virtue of choosing a Gaussian initial distribution, g'^{t) = 
for t < to- This may be less restrictive than it seems since the results are essentially 
/3-independent. 

We now have the explicit form for [$/]:- 



V<l> V(t)^e'^°^^''U r'0+P0_e*(^['^+]-'^['^-l)5[0+(tf) - $fl 

f r'03p0+P0- e^5o[03]+i(s[<^+]-5[<A-]) s[(f)+{tf) - (125) 



where the boundary condition B is 0±(to) = 03(^o) = 03(^o — ^/^o)- This can be written 
as the time ordering of a single field:- 

pt^[^f] = I P0e^^^[<^] 5[0+(t^) - (126) 

along the contour C = C+ ©C_ ©C3, extended to include a third imaginary leg, where 
takes the values 0+, 0_ and 03 on C+, C_ and C3 respectively, for which Sc is >S'[0+], 
5'[0_] and 5*0 [03]. This corresponds to working with a curve of the form shown in figure 
^. We stress again that although 5*0 [0] may look like the quadratic part of S[(j)], its role 
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Figure 4. Closed Time Path with Initial Conditions 

is solely to encode the initial distribution of configurations $ and need have nothing to 
do with the physical action. Henceforth we drop the suffix / on $ j and take the origin 
in time from which the evolution begins as to = 0. 

We perform one final manoeuvre with pt[^] before resorting to further approxima- 
tion. This will enable us to avoid an ill-defined inversion of a two-point function later 



on. Consider the generating functional:- 



Z[j+,j.,js] = j^Vcj) expj^^cM +^ / (127) 
where / is a short notation for:- 

/roo 
3<P^ dt[3+{t)<P+{t)-3_<p.{t)]+ Ut)Mt)dt, (128) 
Jo Jo 

omitting spatial arguments. It may seem from the exponent S'o[03] + — S[(f>^] that 
there is no communication between fields in different branches of C. This is not the case. 
Suppose that ^2 > > ^o- Then the two-point function < T(0_(t2)0+(ii)) > equals 
< 0+(t2)0+(^i)) > since, if we remove the contour for t > t2, ^-(^2) = (p+ih)- Thus, 
in fact , we have a nondiagonal matrix of propagators, of which the 2x2 submatrix 
between the 0+, 0_ fields is the most important. It is for this reason that the propagator 
(3) is incomplete. 

Then introducing aa{x) where a = 1, . . . , A^, we find:- 

Pt^.[^] = J T^aJ^T^<P exp|z5c[0]| expjz J c?^x«,(f)[0+(t;, f) - <l>(f)],|(129) 
= Jva exp|-i J aa<^a^ Z[a, 0,0], (130) 

where a is the source a{t, x) = a{x)6{t — tj). As with V(j), Va denotes Yli Vaa- 

We have seen that analytic progress can only be made insofar as pt[^] is itself 
Gaussian, requiring in turn that 0,0] be Gaussian in the source a. In order to 
treat the fall from the false into the true vacuum, at best this means adopting a self- 
consistent or variational approach i.e. a Hartree approximation or a large N-expansion.. 
However, if we limit ourselves to small times t then Pt[^] is genuinely Gaussian since 
the field has not yet felt the upturn of the potential. That is, we may treat the potential 
as an inverted parabola until the field begins to probe beyond the spinodal point. The 
length of time for which it is a good approximation to ignore the upturn of the potential 
is greatest for weakly coupled theories which, for the sake of calculation we assume, but 
physically, we expect that if the defect counting approximation is going to fail, then it 
will do so in the early part of the fall down the hill. 

7. DEFECT PRODUCTION AT A QUENCH 
7.1. Counting The Defect Density 

To calculate the defect density requires knowledge of pt[^] which, for the moment, we 
write as 

ptm =Mexp(^-^J d'^xd^'y ^aix)Kabix-y;t)^hiy)), (131) 

with Kab = SabK and A/" a normalisation. The circumstances under which a Gaussian 
is valid will be examined later. For weakly coupled theories we shall see that, for short 
times after to at least, a Gaussian pt[$] will occur. If this is taken for granted it has 
been shown by Halperin how to calculate the number density of defects. Postponing 



the calculation of K until then, wc quote those of Halperin's results that are relevant 
to us. 

Suppose that the field (t>{t,x) takes the particular value We noted that 

strings are tubes of false vacuum. Wc count the strings passing through a surface S by 
identifying them with its zeroes in S. Similarly, monopoles are balls of false vacuum in 
which the field vanishes. The only way for a zero to occur with significant probability 
is at the centre of a topological defect so, but for a set of measure zero, all zeroes are 
topological defects. 

Consider the 0{D) theory in D spatial dimensions, with global monopoles. Al- 
though less relevant than strings for the early universe they are slightly easier to perform 
calculations for. Almost everywhere, monopoles occur at the zeroes of labelled 
Xi, i = 1, 2, . . ., at which the orientation $(a;)/|$(a;) | is ill-defined. A topological winding 
number rii — ±1 can be associated with each zero Xi by the rule:- 



Ui = sgndet(9a$6) 



(132) 



Monopoles with higher winding number arc understood as multiple zeroes of $(af) at 
which the rii are additive. The net monopole density is then given by:- 

Pnet{x) ^"^UiSix - Xi). (133) 

i 

The volume integral of this gives the number of monopoles minus the number of anti- 
monopoles. The correlations of pnet give us information on monopole-(anti)monopole 
correlations but, in the first instance, we are interested in the cruder grand totals. The 
quantity of greater relevance to us, is the total monopole density:- 

p(f) = ^5(f-f,), (134) 

i 

whose volume integral gives the total number of monopoles plus antimonopoles. in the 
volume of integration. 

Now consider an ensemble of systems in which the fields $ are distributed accord- 
ing to pt[$] at time t. Then, on average, the total monopole density is:- 



Pm{t) = {P{x))t = <^E^(^-^i)/^' (135) 

where the triangular brackets denote averaging with respect to and the suffix m 
denotes monopoles. That is:- 

= I (136) 

with normalised so that = 1. The translational invariance of the 

Gaussian kernel of the probabililty density ensures that p{t) is translationally invariant. 
In terms of the fields vanishing at Xi, Pm{t) can be re-expressed as:- 

Pm{t) = { S^'lMx)] I det (daM^)) I )t- (137) 



The second term in the brackets is just the Jacobian of the transformation from x to 
$(f). 



It follows from the Gaussian form of the probability density that the are 
individually and independently Gaussian distributed with zero mean, as 

( <^aix) )t = 5abWi\x-y\;t), (138) 

where Ty(|a; — y\;t) = K^^{x — y;t). So also are the first derivatives of the field da^b, 
which are independent of the field:- 

{^c{x)daMy))t = (139) 

due to the fact that W is dependent only on the magnitude oi x — y. 

Thus, the total defect density may be separated into two independent parts:- 

Pmit) = (5^[$(f)])t (|det(9,$6)|)f (140) 

The first factor is easy to calculate, the second less so. 
Consider first the delta-distribution factor:- 

(5[$(fo)])t = /©^^[^(xo)] exp|-i|(i^f(i^y$(f)i^(f-y;t)$(y)| (141) 
= j^ajv^ exp|iQ;$(fo) - 1 / (fx(fy^x)K{x - y; t)^y) |l42) 

where 0{N) indices and integrals over spatial variables have been suppressed and ~ 
da/2Ti. On defining a = S{x — Xo)a, we find:- 

= j ^ajv^ expj j d^x (ia{x)^{x) -\j d^xd^y^x)K{x - y; t)^y)^ | 
= J ^acxpl^-^ J d^xd'^ya{x)W{\x-y\;t)a{y)^ (143) 
= (la expj -i jd^xd^yS{x-Xo)aW{\x-x'\;t)aS{x'-Xo)^ (144) 



27r\y/K^J 27r($$) 27rW{0;t) 



(145) 



Consider now the second factor. Writing out the determinant explicitly for N — 
D — 2, and exploiting the fact that the field is Gaussian, we have:- 



(I det{daM^))\yt 



dct{da(f)b{x)) 



(146) 



^(9i0i9202)' + (dihd2(f>iy - 29i02920i9i0i9202^^.(147) 

The first term may be factorised into a product of two Gaussian variables and calculated 
as follows:- 

((5l0A02)')t = ((5l0l)')t((5202)')t (148) 

= [-Sndid^Wi\x\;t)] [-^s^sl^df |; t)] = [did^W i\x\; t)]%U9) 
where iy(|a;|;t) is as before. Fourier transforming the two-point function, we find:- 

( (9i0i9202)' )t = (d^d, J e'^-^W{k; t) fk^ (^8282 J e'^-^W{k; t) fk^ (150) 

= Q k^cos\e)W{k;t)fk^^ (151) 
= [V^W{0;t)]^ (152) 



A similar result applies for the second term whereas the third term vanishes. The final 
result is:- 

N/2 

(153) 

where the second derivative in the numerator is with respect to x = Cn is l/vr^ 
for N = D = 3 and l/27r for N = D = 2 (had we performed the calculation in D = 2 
dimensions), the difference coming entirely from the determinant factor |p. 

Let us now consider the case of global strings in D = 3 dimensions that arise from 
the 0(2) theory. Strings are identified with lines of zeroes of 0()f:, x) = $(x) and the net 
vortex density (vortices minus antivortices) in a plane perpendicular to the z-direction 
is:- 

Pnet,i(x) = S'l^x)] e,,kid,^,)idk^2), (154) 

with obvious generalisations to = -D — 1, for all N in terms of the Levi-Cevita symbol 
^ii,i2,...,iD- As before, the total vortex density is of more immediate use. On a surface 
perpendicular to the i-direction this is:- 

p,(f) = S^mx)] \e,,k{dj^i)idk^2)\. (155) 

in analogy with the monopole case. The expectation value of this total density, when 
calculated as before reproduces the same expression (where s denotes string s):- 

N/2 

(156) 

independent of i, but for N = D — 1. Thus, whether we are concerned about global 
monopole or global string density, once we have calculated iy(|x|;t) we can find the 
total string density. Similar results apply for the correlations between net densities, 
which are important in determining the subsequent evolution of the defect network. 

7.2. Evolution Of The Defect Density 

The onset of the phase transition at time t = is characterised by the instabilities 
of long wavelength fiuctuations permitting the growth of correlations. Although the 
initial value of (0) over any volume is zero, the resulting phase separation or spinodal 
decomposition will lead to domains of constant (0) whose boundaries will evolve so 
that ultimately, the average value of in some finite volume, will be non-zero. That 
is, the relativistic system has a non-conserved order parameter. In this sense, the 
model considered here is similar to those describing the transition in liquid helium 
or transitions in a superconductors. 

Consider small amplitude fiuctuations of 0a, at the top of the parabolic potential 
hill described by V{^) = ^m'^{t)(j)l. At t < 0, m^{t) > and, for t > 0, m'^{t) < 0. 
However, by t ~ At, m?{t) and g'^ have achieved their final values, namely — /i^ and g'^. 
Long wavelength fiuctuations, for which < —rn^it), begin to grow exponentially. If 

their growth rate \J —7rfi{t) — [fc^ is much slower than the rate of change of the 

environment which is causing the quench, then those long wavelength modes are unable 
to track the quench. For the case in point, this requires mAt -C 1. We take this to be 
the case. To exemplify the growth of domains and the attendant dispersal of defects. 



pm{t) — Cn 



iy"(0; t) 



W{Q-t) ' 



Ps,i{^) =< pi{x)>= C, 



W"{0;t) 
W{0;t) 



it is sufficient to take the idealised case, At = in which the change of parameters at 
t = is instantaneous. That is, m^(t) satisfies 



> 
-/i2<0 



if t < 0, 
if t > 



(157) 



where for t < 0, the field is in thermal equilibrium at inverse temperature Pq. As for 
g'^it), for t < 0, it has already been set to zero, so that Ptol^] be Gaussian. For small 
t, when the amplitude of the field fiuctuations is small, the field has yet to experience 
the upturn of the potential and we can take g'^{t) = then as well. At best, this 



can be valid until the exponential growth 



/ie'^* in the amplitude reaches the 



point of infiection |0| ^ fi/y/g^, that is fit ~ 0(ln(l/(yf^)). The smaller the coupling 
then, the longer this approximation is valid. As noted earlier, it should be possible 
to perform more sophisticated calculations with the aim of evolving the defect density 
right through the transition. For our present purposes, however, the small time or 
Gaussian approximation is adequate. 

We are now in a position to evaluate identify K and calculate the defect 

density accordingly. Sc[<p] becomes S'o[03] on segment C3 so setting the boundary 
condition 0+(O,x) = 03(0, x) = (j)s{—iPo,x) and we have 



d X 



158) 



on C+. The Gaussian integrals can now be performed to give 



= J'^Ci exp|-z J rf^a;aa$a|exp|^y" £x £y aaix)G{x - y;t,t)abiy)y (159) 

where G{x — y;t,t) is the equal time correlation, or Wightman, function with ther- 
mal boundary conditions. Because of the time evolution there is no time translation 
invariance in the double time label. As this is not simply invertible, we leave the a 
integration unperformed. The form is then a mnemonic reminding us that = G. 

In fact, there is no need to integrate the as since from the previous equation it 
follows that the characteristic functional (exp{z / Ja^a})t is directly calculable as 



expjz j Ja^aY)^ = Jl^^Ptm expji J Ja'^a 



= ^^^[2] d^y ja{x)G{x ~ y]t,t)jai 
Thus for example, the first factor in the monopole density pm(t) is 



(160) 
(161) 



J dj expi^^jlG{0;t,t) 



-iG{0;t,t)] 



-D/2 



(162) 
(163) 



with suitable normalisation, without having to invert G{0;t). Thus, on identifying 
—iG{x;t,t) with W{x,t) as defined earlier, Pm(t) becomes H 



Pmit) = C, 



N 



-iG"{0;t,t) 
-iG{6;t,t) 



N/2 



(164) 



where —iG{x; t, t) has to be calculated from the equations of motion, subject to the 
initial condtition. 

Details are given by Boyanovsky et al., and we quote their results, which give 
—iG{x;t,t) as the real, positive quantity 



—iG{x; t, t) 



ik.x 



2cj<(A;) 
l + Ak{cosh{2W{k)t) - 1 



coth(/5ou;<(A;)/2)x 



(165) 



1 + afc(cos(2u;>(A;)t) — 1 



with 



Wl{k) 



Ifcp + ml 



.2 



A, 



- 1+^ , 

2 I W\k) / 



(166) 
(167) 
(168) 

(169) 



Oik 




(170) 



The first term is the contribution of the unstable long wavelength modes, which relax 
most quickly; the second is that of the short wavelength stable modes which provide the 
noise. The first term will dominate for large times and even though the approximation 
is only valid for small times, there is a regime, for small couplings, in which t is large 
enough for cosh(2/it) ?s |ea;p(2/it) and yet /it is still smaller than the time 0(ln 1/ g^) at 
which the fluctuations sample the deviation from a parabolic hill. In these circumstances 
the integral at time t is dominated by a peak in the integrand /j-D-ig2vy(fc)t around 



kc, where 



tkl 



[D-l] 



/i l + 0( 



(/it))' 



(171) 



The effect of changing /?o is only visible in the 0(1/ jjt) term. In the region |x| < ^Jt/n 
the integral is dominated by the saddle-point at kc, to give 



iG{x;t,t) = W{x;t) ^ W{0;t) expl 



-fix 



X 



smc 



t//i. 



for D = 3, where 



W{0]t) ^ c- 



(172) 



(173) 



(/it)3/2^ 

for some C, which we don't need to know. The exponential growth of G{0] t) in t reflects 
the way the field amplitudes fall off the hill ($) = 0. It is sufficient for our purposes to 
retain D = ?> only. 

After symmetry breaking to 0{N — 1) the mass of the Higgs is mn = v^/t with 



cold correlation length ^(0) = . On identifying e 



■fj.x^/8t 



as e 



we interpret 



at) = (8v^)i/2/u(o)' 



(174) 



as the size of Higgs field domains. This t^/^ growth behaviour at early times is char- 
acteristic of relativistic systems (with a double time derivative) with a non-conserved 
order parameter. 




To calculate the number density of defects at early times we have to insert this 
expression for —iG or W into the equations derived earlier. Expanding W{x; t) as 

*) = H/(0; t) exp(^) ( 1 - 1^ + o(^) ). (175) 

gives [H 

for an 0(3) theory with monopoles in three dimensions and 

(177) 

for an 0(2) theory with strings in three dimensions. The first observation is that the 
dependence of the density on time t is only through the correlation length ^(t). As the 
domains of coherent field form and expand, the interdefect distance grows accordingly. 
This we would interpret as the domains carrying the defects along with them on their 
boundaries. Secondly, as we signalled in the introduction, there is roughly one defect 
per coherence size, a long held belief for whatever mechanism. However, in this case 
the density is exactly calculable. It is also possible to use Halperin's results to calculate 
defect-defect correlation functions. Details will be given elsewhere. 

Finally, we mentioned in our introduction that superfiiud He is, if anything, more 
likely to provide an environment in which these mechanisms can be tested. One way to 
convert the relativistic formalism above into that for a non-relativistic field theory is to 
include a chemical potential of value m equal to the rest mass of the quanta. Provided 
f3m 3> 1 antiparticles are eliminated and effective particle energy is non-relativistic. 
Otherwise everything goes through as before. Results will be presented elsewhere. 



8. LATER EVOLUTION 

Once field fiuctuations are large enough to reach the points of inflexion of the classical 
potential V{(j)) = g'^{<f)'^ — rj^Y/S, the simple Gaussian approximation breaks down, if 
not before. For a while the evolution of domains is controlled by the dissipation of the 
field energy in particle production. However, because of their topological stability some 
defects will survive. 

One way to proceed is to try to develop Landau-like Langevin equations for the 
evolution of coarse-grained fields in V{(f)). From our previous considerations the relevant 
scale would be ^(t,), where tj = 0{ln{l/ g^)) is the time at which the field fiuctuations 
begin to feel the points of infiexion of V^(0). This is essentially the approach of Davis 
and Martin |37], although they choose ^g-, the Ginzburg correlation length as their 



starting point. From the point of view of the equations this does not matter, and this 
is all that concerns us here. As yet we have no new results to report, but we conclude 
these lectures with a few comments on dissipation that might help in developing a 
programme along the lines of |]3^. In large part, we are just paraphrasing the most 



recent work of Boyanovsky et al. [§]. 

It is a truism to say that coarse-graining introduces dissipation into the otherwise 
time-reversal invariant evolution equation for the fields. To see how this comes about, 
we begin with the simplest example of coarse-graining, the effective potential. 



8.1. The Effective Potential 

We have discussed the effective potential earlier since, as the energy density, it deter- 
mines the ground state of the system, and hence its phase. However, we avoided a 
detailed calculation of it. For exemplary purposes it is sufficient to consider a single 
scalar field at zero temperature. 

The effective potential is obtained most succinctly by summing over field 

histories in the partition function Z, subject to the constraint 

d'x4>{t,x) = 4>. (178) 

where Q is the large space-time volume of the system. This total coarse-graining of 0, 
both in space and time, precludes any evolution equation, leading just to 



d(f) 



(0) = (179) 



as the definition of the ground state of the theory. Nonetheless, eq.(179) is the simple 
counterpart of the Langevin equation that we shall attempt later. 

To see how (179) occurs it is most useful to make the separation 

(p{t,x) = ^ + r){t,x), (180) 

where rj contains no fiuctuations with k — 0. That is, 

J d'^xr]{t,x) = 0. (181) 

To calculate we begin with the generating functional Z\j], in which is coupled 
to a constant source j as 

Z\j] ^ jT>(f) e^^[<^]+^j /^ (182) 

where is given in (4). For the sake of argument is taken positive. Upon the 
decomposition (180) Z[j] becomes 



Z\j] ^ J d^ J Vt] S i^J Tjj e*^[^+'']+^^'"'^. (183) 
A series expansion of S[(j) + rj] in powers of r] leads to 

Z[j] = J d4>e'^^^-'^^^^^ J'^^^iJ e'^^'^'"'', (184) 

where 

SefAv] = I d'x {^d^vd^v - IK + sg'^V'^W + •••) • (185) 

The constraint (181) eliminates all terms in Seff linear in rj. To one-loop it is only 

necessary to 

i. retain quadratic terms in rj in Sgff 

ii. ignore the 5-function. 



On performing the 7] integration we find 
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Z[j] - J e'^^^'^-^^ff^'^^\ (186) 

where 

Veffi^) = V{^) + \i j d^k ln{-P + m^ + 3^202/2). (i87) 

is the effective potentiaL In the Q ^ oo hmit the integral is dominated by a stationary 
phase for satisfying (179), which thereby determines the ground state of the theory. 

We conclude this introductory section by making some simple observations on 
approximations. The one-loop approximation is correct to 0{h) and all orders in g^. 
Its 5f^-expansion, to second order, is 

K//(0) = V{<P) + 1 S + 1 / + Oig% (188) 

where E is the one-loop self-mass and I is the one-loop two-vertex function. 

If we had only been given (188) as a first approximation to Veff we would know 
that it would need to be elevated to the one-loop potential (187). Further, in many 
circumstances even this is not sufficient. Any attempt to incorporate non-perturbative 
effects in V^ff would require that, at least, we extend the one-loop potential (187) 
to a self-consistent one-loop potential, using the Hartree approximation or a large-N 
calculation of the 0{N) theory. For the case in hand, when N is directly related to 
spatial dimension D as N = D or N = D — 1 the latter is not useful. 

When we attempt a Langevin equation it is the counterpart of (188) that we can 
most simply derive. We shall do this in the knowledge that both the steps listed above 
need to be implemented before we can have believable results. 



8.2. Spatial Averaging 

As the next step we determine the evolution equation for the spatial average 

<l>(t) = lim 0Jt) (189) 

for (l)v(t) of (43). Let us return to the closed time-path of (127) where, as there, we have 
assumed an initial thermal distribution of states at inverse temperature /Sq. As with 
the effective potential we wish to separate the field into those modes that describe 
the 'system', and those that describe the 'environment' that will be integrated over. In 
this case these latter are the k modes, permitting the decomposition 

(j)±{t,x) = <^±{t)+r]±{t,x), (190) 

where 



d^xr]±{t,x) = (191) 

at all times t. Although $+(t) ^ $-(i), < >=< $-(t) > since each is equally 

acceptable as the physical field. They only differ by the A; = fluctuations, 

a(t) = $+(t)-$_(t). (192) 

The effective coarse-grained field is 

<l>(t) = i($+(t) + $_(t)) (193) 



with common expectation value. 

The generating functional for the effective field $ and its correlations is 



Z[j]= r'03r'0+r'0_e*^«['^^l+^(^['^+]-^t'^-l)e^''/^'*, (194) 
J c 

where 5*0 sets the initial conditions and, from above 

0±(t, x) = <^{t) ± ia(t) + r/±(t, x). (195) 

On substituting in Z, the constraint (191) guarantees that there are no linear terms 
in T]±, just as no terms in rj survived in (185). Thereafter, at one loop we can ignore 
the constraint in the path integral. Further, to lowest nontrivial order it is sufficient to 
neglect terms 0{ri^) in 7]^ and terms O(cr^) in a. As a result 

S[<^ + + r]+] - S[<l> - la + T]^] = (196) 

/ dtdx (a(t)^ + [^idr^^f - iMlr^l] - [^idr^.f - iM!^]) (197) 

where 

Ml=m^ + ^/($ ± |(t)2 (198) 
and 1. 1 1 

£($) = _ ^^2^2 _ ^^2^4_ ^^gg) 

2 2 8 

We can now integrate out the ?7±, remembering that they are coupled and expand the 
resultant trace of the logarithm (cf. (187)) in powers of g^. At second order in the 
effect is to give Z[j] as 

Z[j] = J V^Va e*-[^«//[*''^l+/^*l, (200) 

where Sejf is 

= - Jdt (^it)L{^{t)) + \i j dtdt' a{t)K{t, t'; $)a(t') (201) 
for L and K real. L{^) is the dissipative form 



(202) 



the counterpart to dVejf/dcj) for Veff of (188), in which S(t) is the 'tab' diagram coun- 
terpart to S .I(t,t'), real anc? retarded, is built from one-loop two- vertex diagrams, the 
counterpart to / of (188). With propagators G±± for r]±, more one-loop diagrams are 
possible than for the effective potential. Details are given in Boyanovsky et al. but 
see also [|39|. At one-loop I(t,t') provides dissipation. Meanwhile K{t,t'), of the form 
{3g'^^{t)/2)R{t,t'){3g'^^{t')/2), is also real, where R is the other combination of one- 
loop two- vertex diagrams, but gives a relative imaginary part to Sefj- K describes the 
'noise' of the environment, most easily seen by using the identity 

e-bS'^K'^ = j Vi 6^5^/^^''^+™/^'' (203) 



to define P[^, $] by 

P[^,$] = e-^^/«^"'«. (204) 

Then Z[j] becomes 

Z[j] = J V^V^Va P% $]e-^''/'^(^-^)e'^/^*. (205) 

The (T-integration is trivial, giving 

Z[3\ = J V^Vi P% $]5[L($) - e]e^''/^■^ (206) 

where the square brackets denote a (5-functional, valid at each time. That is, for each 
noise function ^(t), $(t) satisfies the dissipative Langevin equation 

L{m) = m (207) 

for L of (201). Both S(t) and I{t,t') depend on the initial conditions and describe 
the shift in mass and the dissipative effect of the coarsegraining from a particular 
distribution of initial states. As it stands, if the natural frequencies (masses) of the 
initial distribution do not match the frequencies (masses) of the potential there will be 
transient shocks. However, since the expression (201) can only be valid for short times, 
at best, 'transient' effects may be important. The noise has distribution P[^,$[^]], 
depending on through $ but, in the Gaussian approximation, $(t) =< $(t) > satisfies 
the purely dissipative equation 

L($(t)) = 0. (208) 

Equation (207), an Ehrenfest equation for field theory, can be obtained directly form 
the functional Schroedinger equation, if required. 

As we noted with the effective potential, this approximation has to be elevated 
to a full one-loop equation and then to a consistent one-loop approximation before it 
begins to be believable. Nonetheless, this shows how dissipation and noise arise in 
quantum field theory. [It should be remembered that, in an expanding universe, the 
evolution of the metric provides its own dissipation, 3H^]. 

As yet we are not in a position to build upon the understanding of domain for- 
mation that we have developed in the previous sections. Complete spatial averaging 
removes any reference to defects in the same way that the space-time averaging of the 
effective potential obliterates any information of finite-size fluctuations. [This can be 
seen very easily in the context of equilibrium theory, where the effective potential es- 
sentially describes the (logarithm of) the probability that a particular field average is 
achieved over all space-time [ED| . 

In order to describe defects it is necessary to coarse-grain on correlation-length 
scales. That is, we separate fields (p{t,x) into 'system' and 'environment' as 

0(t, x) = (pfit, x) + 0i_/(t, f), (209) 

where 

^f{t,x) = I cffc3e^^^-%(t)/(fc) (210) 
contains only wavelengths with with < for some ^. That is. 



fik) = e{h-\k\), 



(211) 



where kc = We can now try to integrate out the short-wavelength modes, in 

the first instance in a one-loop approximation at There will be dissipation and 

noise as before, arising from diagrams with similar ultraviolet behaviour, but different 
infra-red properties. The resulting Langevin equation should permit us to follow the 
evolution of defects, but serious sums have yet to be done. This seems a good place to 
stop. 

9. CONCLUSIONS 

In these notes we have made several attempts to understand defects produced at a phase 
transition. We will only recapitulate our main result. This is that defect density at 
formation is controlled by the field correlation lengths essentially in the way anticipated 
by Kibble and others. 

It is one thing to count defects in fluctuations, another to demonstrate how they 
freeze out as the system becomes cold. For quasi-equilibrium in the Ginzburg regime 
defects appear copiously, but we are unable to show how they can survive as temperature 
is lowered. More hopefully, defects produced at a quench by phase separation do track 
domain boundaries and have every likelihood of freezing-in in numbers determined by 
these domains, but as yet we do not have the machinery to follow them through. Much 
more work remains to be done. 
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